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Methanol  autothermal  reforming  was  thermodynamically  analyzed  using  FLUENT  software.  The  calcu¬ 
lation  methodology  using  this  software  is  simple  and  convenient,  and  its  validity  was  confirmed  by 
comparing  the  obtained  data  with  previous  studies.  As  a  function  of  the  effects  of  temperature,  pressure, 
molar  steam-to-carbon  ratio  (S/C),  and  molar  oxygen-to-carbon  ratio  (O/C)  on  the  objective  products, 
favorable  operational  parameters  were  evaluated,  under  which  H2  yield  maximizes,  the  CO  molar  fraction 
minimizes  and  carbon  deposition  can  be  eliminated.  The  equilibrium  constants  of  the  possible  reactions 
involved  in  oxidative  methanol  steam  reforming,  coupled  with  the  reaction  mechanism  for  the  entire 
investigated  temperature  range,  were  elucidated  and  discussed.  On  the  basis  of  the  concluded  possible 
mechanisms,  three  areas  are  inferred.  In  each  individual  area,  H2  or  CO  yield  reached  a  maximum,  or 
solid  C  was  efficiently  suppressed.  Therein,  a  favorable  operational  range  is  proposed  to  assure  the  most 
optimized  product  yield. 

©  2008  Elsevier  B.V.  All  rights  reserved. 


1.  Introduction 

Increasing  environmental  awareness  and  the  associated  devel¬ 
opment  of  fuel  cell  technologies  has  directed  greater  attention  to 
fuel  processing  technologies  for  hydrogen  generation.  Presently, 
hydrogen  refueling  stations  are  rare,  and  do  not  have  the  ubiqui¬ 
tous  infrastructure  available  to  gasoline,  diesel-fuel,  and  natural  gas 
fueling  technologies.  Furthermore,  the  complex  compositions  of 
these  petroleum  products,  specifically  their  sulfur  content,  hinder 
their  application  as  feedstocks  for  generating  hydrogen,  especially 
in  countries  like  China.  Natural  gas  has  a  higher  hydrogen-to- 
carbon  ratio  (H/C),  which  is  more  favorable  for  producing  hydrogen; 
however  it  is  only  used  in  on-site  hydrogen  generation  due  to  its 
inconvenient  storage  and  transportation. 

Significant  effort  has  been  devoted  to  the  development  of  liquid 
fuels  for  on-board  hydrogen  generation,  especially  liquids  com¬ 
posed  of  hydroxyl-containing  hydrocarbons.  Methanol  in  particular 
is  regarded  as  the  most  feasible  fuel  for  on-board  hydrogen  genera¬ 
tion,  primarily  due  to  its  higher  H/C  and  low  reforming  temperature. 

Three  technologies  have  historically  been  used  to  optimally  gen¬ 
erate  hydrogen  from  methanol,  namely,  partial  oxidation  (POX), 
steam  reforming  (SR),  and  oxidative  steam  reforming  (OSR)  or 
autothermal  reforming  (ATR).  POX  is  an  exothermic  process  that 
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requires  a  high  reaction  temperature,  and  unfortunately  produces 
products  with  low  hydrogen  concentrations.  SR  is  an  endother¬ 
mic  process,  and  requires  additional  heaters  or  burners  to  provide 
driving  energy  for  the  reaction.  Exothermic  processes  have  an 
intrinsically  faster  transient  response  compared  to  endothermic 
processes.  It  is  important  to  note  that  ATR  combines  the  partial 
oxidation  and  steam  reforming  processes,  wherein,  by  properly 
adjusting  the  reactive  conditions,  the  reaction  system  can  be  run 
autothermally.  The  reaction  temperature  and  the  quality  of  the  fuel 
cell  feed  in  the  ATR  lie  between  POX  and  SR.  Further,  a  higher  ther¬ 
mal  efficiency  is  possible  with  an  optimally  developed  ATR  process. 

For  such  a  complex  reaction  system,  thermodynamic  equilib¬ 
rium  calculations  can  be  performed  via  two  common  approaches 
[1  ].  One  approach  is  based  on  chemical  equilibrium  constants,  and 
is  called  the  chemical  stoichiometric  method,  while  the  other  is 
derived  from  the  minimization  of  the  Gibbs  free  energy.  The  for¬ 
mer  is  inapplicable  to  liquid  or  solid  species  in  the  reaction  system. 
For  a  heterogeneous  system,  the  approach  using  the  minimization 
of  the  Gibbs  free  energy  is  commonly  used  for  thermodynamic 
analysis.  Commercial  software,  based  on  the  Gibbs  free  energy 
function,  has  been  widely  used  for  thermodynamic  calculations, 
such  as  AspenPlus™  [2]  and  HSC-Chemistry  (Outokumpu,  Finland) 
[3].  Equilibrium  calculations  can  also  be  performed  through  pro¬ 
gramming. 

Chan  [4]  thermodynamically  analyzed  ATR  based  on  the  For¬ 
tran  code  developed  by  NASA  under  the  assumption  that  the  inlet 
temperature  of  the  reactants  was  298.15  K  and  the  fuels  were 
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completely  converted.  Generally,  this  scenario  is  not  true  in  most 
practical  cases,  therefore  we  used  optimized  practical  operational 
parameters  dependant  on  a  similar  method  described  in  [5];  how¬ 
ever  the  solutions  and  the  data  post-processing  are  complicated. 
Win  et  al.  [6]  investigated  a  complex  numerical  algorithm  for 
methanol  steam  reforming.  The  probability  of  carbon  formation 
in  practical  scenarios  was  determined  by  the  direct  inclusion  of 
carbon  in  the  objective  function  of  the  minimization  scheme,  such 
that  an  intricate  solution  for  the  minimization  of  the  Gibbs  free 
energy  function  was  unavoidable.  The  stoichiometric  method  was 
also  used  to  perform  thermodynamic  analysis  [7];  however,  the 
assumption  of  no  formed  liquid  or  solid  species  limits  its  appli¬ 
cation. 

In  this  work,  a  method  has  been  developed  to  perform  the 
thermodynamic  analysis  of  methanol  OSR  for  PEMFC.  Independent 
parameters  that  effect  ATR  reactor  performance,  that  is,  the  inlet 
feed  temperature,  pressure,  S/C,  and  O/C  were  investigated.  As  a 
function  of  the  equilibrium  constants  of  the  possible  reactions  in 
the  reactive  system,  the  formation  mechanisms  of  the  object  prod¬ 
ucts  were  established.  On  the  basis  of  the  thermodynamic  analysis 
and  formation  mechanism,  the  optimal  operational  parameters 
that  maximized  H2  yield,  minimized  the  CO  molar  fraction,  and 
eliminated  carbon  deposition  were  confirmed.  In  addition,  the  pro¬ 
posed  analysis  method  was  validated  by  comparing  the  obtained 
results  with  those  reported  in  the  literature. 


2.  Theorem  and  simulation  methodology 

2.1.  Theorem 


In  contrast  to  the  chemical  stoichiometric  method,  the  mini¬ 
mization  of  Gibbs  free  energy  method  is  applicable  to  systems  with 
condensed  phase  species.  H2  generation  systems  using  a  fossil  fuel 
feed  can  produce  coke  under  certain  thermodynamic  conditions. 
To  avoid  the  negative  impact  of  coke  formation,  these  thermody¬ 
namic  conditions  can  be  identified  using  the  minimization  of  Gibbs 
free  energy  method.  A  closed  system  at  equilibrium  satisfies  the 
minimized  Gibbs  free  energy  for  a  specified  temperature  and  pres¬ 
sure,  while  is  also  simultaneously  constrained  by  the  conservation 
of  atoms.  In  this  vein,  many  equations  have  been  established  and 
solved  using  species  composition  and  thermodynamic  data.  A  sig¬ 
nificant  amount  of  code  is  available  that  can  solve  these  equations, 
most  notably  the  NASA  and  STANJAN  equilibrium  codes.  NASA  pro¬ 
gram  has  been  used  in  the  analysis  of  generating  hydrogen  process 
fed  by  fossil  fuels  [4,8].  In  addition,  STANJAN  has  been  adopted 
in  some  commercial  thermodynamic  software,  such  as  CHEMKIN. 
The  Gibbs  free  energy  minimization  method  and  the  STANJAN  code 
were  applied  to  preprocess  the  non-premixed  combustion  model 
in  the  FLUENT  software.  A  conserved  scalar  quantity,  known  as  the 
mixture  fraction/,  was  introduced  to  perform  thermodynamic  equi¬ 
librium  analysis  in  the  solving  process.  The  mixture  fraction  can  be 
written  in  terms  of  the  atomic  mass  fraction  as 


f  = 


Zj  -Zj,ox 
Z[fuel  ~  Z[  qx 


(1) 


where  Zx  is  the  elemental  mass  fraction  for  some  element  i.  The 
subscript  “OX”  denotes  the  value  at  the  oxidizer  stream,  and  the 
subscript  “fuel”  denotes  the  value  at  the  fuel  stream.  Due  to  unique 
relations  between  all  of  the  thermochemical  scalars  (species  mole 
fraction yit  density,  and  temperature)  and  the  mixture  fraction,/, 
the  species  mole  fraction  can  be  solved  based  on  the  minimiza¬ 
tion  of  Gibbs  free  energy  using  the  STANJAN  code  [9].  The  species 
mole  fractions  were  then  converted  into  a  dry  basis  and  displayed 


in  figures  by  the  expression: 


y  i,dry 


y\ 

i  -yH2o 


(2) 


In  addition,  the  mixture  fraction,  /  was  transferred  into  the  hori¬ 
zontal  coordinate  of  the  O/C  in  terms  of  the  expression: 


O 
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0.21  — jrr^— 
MAir 


Mq2  \ 
Mch3oh  J 


(3) 


where  MCh3oh>  MAir,  and  Mq2  are  the  CH3OH,  air,  and  02  molecular 
weights,  respectively. 


2.2.  Thermodynamically  feasible  reactions  and  products 

On  the  basis  of  thermodynamic  analysis  and  experimental 
observation  [10-17],  the  following  reactions  are  primarily  involved 
by  neglecting  the  intermediates  and  products  from  side  reactions; 
however,  the  contribution  to  the  product  distribution  by  each  reac¬ 
tion  differs  depending  on  the  specified  operational  conditions. 
Among  them,  the  primary  reactions  contributing  to  the  conversion 
of  methanol  include: 


Methanol  decomposition: 

CH3OH(g)  =  CO  +  2H2,  AH298=  90.67  kj  mol”1  (1) 

Methanol  partial  oxidation  (POX): 

CH3OH(g)  +  0.502(g)  =  C02  +2H2, 

AH298  =  -192.3  kj  mol"1  (2) 

Methanol  combustion: 

CH3OH(g)  +  1.502=  C02  +  2H20(g), 

AH298  =  -675.97  kj  mol”1  (3) 

Methanol  steam  reforming  (SR): 

CH3OH(g)  +  H20(g)  =  C02  +3H2,  AH298  =  49.5  kj  mol"1  (4) 

The  water-gas  shift  reaction  (WGS)  is  favorable  for  generating 
further  H2  at  the  expense  of  CO,  and  occurs  as 

CO  +  H20(g)  =  C02  +  H2,  A H298  =  -41  kj  mol”1  (5) 

When  02  is  in  excess,  the  CO  will  be  transformed  into  C02,  and  a 
partial  depletion  of  H2  will  take  place  via  the  preferential  oxidation 
reaction: 

CO  +  0.502  =  C02,  A H298=  -283  kj  mol"1  (6) 

H2  +  0.502=  H20(g),  AH298=  -241.8  kj  mol"1  (7) 

Solid  C  formation  can  also  occur  via  the  following  reactions  [10]: 

2CO  =  C(s)  +  C02,  AH298=  -  172kJmor1  (8) 

CH4=  C(s)  +  2H2,  AH298  =  75.6  kj  mol”1  (9) 

CO  +  H2=  C(s)  +  H20(g),  AH298=  -131.3  kj  mol"1  (10) 

C02  +2H2=  C(s)  +  2H20(g),  AH298=  -90.13  kj  mol"1  (11) 

If  the  molar  S/C  is  kept  lower,  methanation  reactions  can  also  occur: 

CO  +  3H2  =  CH4+H20(g),  AH298=  -206.15  kj  mol”1  (12) 

C02+4H2=  CH4  +  2H20(g),  AH298  =  -leskjmor1  (13) 

In  addition,  it  is  thermodynamically  possible  to  form  ammo¬ 
nia  and  hydrogen  cyanide  in  the  presence  of  N2;  however  their 
respective  effluent  fractional  concentrations  are  exceedingly  small 
(in  the  range  of  ppm  or  ppb)  under  the  simulation  conditions 
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[11],  and  are  therefore  neglected.  As  per  the  analysis  above,  H2, 
CO,  C02,  N2,  H20,  02,  CH3OH,  CH4,  and  solid  graphite  carbon 
C(s)  are  all  the  plausible  products  resulting  from  the  OSR  of 
methanol.  The  impact  of  various  operational  parameters  on  these 
equilibrium  compositions  and  the  adiabatic  temperature  were  per¬ 
formed  using  this  methodology.  In  addition,  the  validation  of  the 
methodology  was  also  confirmed  on  the  basis  of  the  referred 
results. 

3.  Results  and  discussion 

3.1.  Validation  of  the  methodology 

To  confirm  the  reliability  of  the  methodology,  the  FORTRAN  code 
developed  at  the  NASA  Lewis  Research  Center  was  also  used  to 
iteratively  compute  the  chemical  equilibrium  products.  The  FOR¬ 
TRAN  program  is  also  based  on  the  minimization  of  Gibbs  free 
energy,  but  uses  a  different  thermodynamic  database  than  the  FLU¬ 
ENT  software.  The  FORTRAN  code  was  compiled  using  the  NASA 
data  format;  however,  The  CHEMKIN  polynomial  format  was  used 
in  the  FLUENT  software.  Thermodynamic  equilibrium  calculations 
were  performed  using  the  two  methods,  respectively.  Equilib¬ 
rium  compositions  and  temperatures  are  showed  in  Figs.  1-5, 
wherein  A  denotes  the  results  of  the  discussed  method,  and  B, 
the  results  from  the  literature.  Partial  results  based  on  the  FOR¬ 
TRAN  code  were  published  in  [5].  From  Fig.  1,  it  can  be  observed 


Fig.  1.  Effect  of  the  O/C  and  S/C  on  the  adiabatic  temperature  (inlet  temperature  of 
the  reactants  was  573.15 1<  at  a  pressure  of  0.1  MPa). 


Fig.  2.  H2  mole  fraction  as  a  function  of  the  O/C  and  S/C  (inlet  temperature  of  the 
reactants  was  573.15  K  at  a  pressure  of  0.1  MPa). 

that  similar  data  trending  validates  the  feasibility  of  the  pro¬ 
posed  methodology.  In  addition,  the  thermodynamic  analysis  of 
the  methane  reforming  process  provides  the  evidence  for  its  appli¬ 
cability  [12].  A  distinct  divergence  occurred  for  coke  formation 
prediction,  which  is  due  to  the  discrepancies  between  the  ther¬ 
modynamic  parameters  supplied  for  solid  C  from  the  two  different 
databases. 

3.2.  Influence  of  inlet  temperature  on  the  equilibrium 
compositions  and  adiabatic  temperature 

The  profiles  of  the  equilibrium  compositions  and  temperature 
were  plotted  as  a  function  of  the  O/C  and  inlet  temperature  in  Fig.  6. 
These  profiles  suggest  that  a  higher  inlet  temperature  favors  FI2 
generation  and  C  removal;  however,  higher  inlet  temperature  also 
affords  more  favorable  conditions  for  CO  formation.  In  the  vein  of 
thermodynamic  analysis,  the  equilibrium  state  generally  depends 
on  the  equilibrium  temperature  of  the  autothermal  system.  There¬ 
fore,  diverse  inlet  temperatures  will  shift  the  equilibrium  state  due 
to  their  indirect  effects  on  the  equilibrium  temperature,  as  can  be 
observed  in  Fig.  6.  At  1000  K,  FI2  yield  maximizes,  and  the  corre¬ 
sponding  O/C  is  between  0.21  and  0.37  for  the  investigated  inlet 
temperature.  Similarly,  the  maximum  CO  molar  fraction  and  C(s) 
removal  respectively  match  an  equilibrium  temperature  of  about 
1300  K  and  1200  K  when  the  S/C  is  1.5.  Operational  pressure  and  the 
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O/C  can  also  influence  the  equilibrium  compositions,  and  are  the 
reasons  that  the  equilibrium  state  is  a  function  of  the  temperature, 
pressure,  O/C,  and  S/C. 


3.3.  Influence  of  operational  pressure  on  the  equilibrium 
compositions 

When  there  are  gases  in  the  reaction  as  reactants  or  prod¬ 
ucts,  a  change  in  the  pressure  can  shift  the  reaction  direction 
due  to  non-conservation  of  chemical  stoichiometric  coefficients. 
For  H2  generation,  many  gaseous  reactions  are  involved  in  the 
system,  and  therefore,  it  is  necessary  to  optimize  the  opera¬ 
tion  pressure.  The  impact  of  different  operation  pressures  on 
the  main  parameters  was  simulated,  as  depicted  in  Fig.  7.  The 
simulation  results  indicate  that  equilibrium  compositions  were 
primarily  constant  when  the  O/C  was  more  than  0.32  at  dif¬ 
ferent  operational  pressures.  When  the  O/C  was  less  than  0.32, 
the  H2  molar  fraction  decreased  slightly  with  increasing  pres¬ 
sure,  and  trended  conversely  for  CO.  Meanwhile,  higher  pressures 
were  favorable  for  C(s)  removal.  In  consideration  of  safety  and 
economy,  the  system  should  be  operated  at  atmospheric  pres¬ 
sure. 


Fig.  3.  Effect  of  different  O/C  and  S/C  on  generated  hydrogen  yield  (preheated  tem¬ 
perature  of  reactants  was  573.15 1<  at  a  pressure  of  0.1  MPa). 
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Fig.  4.  Influence  of  the  O/C  and  S/C  on  the  mole  fractions  of  CO  in  the  reformer  (inlet 
temperature  of  the  reactants  was  573.15  K  at  a  pressure  of  0.1  MPa). 


3.4.  The  optimal  operation  parameter  for  OSR  of  methanol 

When  generating  H2  for  PEM  fuel  cell  applications,  the  pri¬ 
mary  goal  is  to  produce  the  largest  amount  of  FI2  with  a  minimal 
amount  of  CO  effluent  with  a  molar  fraction  in  the  ppm  range. 
To  avoid  catalyst  deactivation  from  C  deposition,  coke  formation 
needs  to  be  suppressed.  Therefore,  the  product  distributions  are 
presented  under  specified  thermodynamic  processing  conditions, 
and  are  aimed  to  maximize  FI2  yield,  while  minimizing  CO  and  coke 
production. 

By  optimizing  the  S/C  and  O/C,  it  is  possible  to  shift  the  reaction 
route  and  change  the  ratios  of  the  possible  occurring  reactions.  As  a 
result,  the  equilibrium  composition  and  temperature  will  approach 
the  expectation  about  the  maximum  FI2  yield,  minimum  CO  and 
coke  formation.  Therefore,  the  impact  of  different  S/C  and  O/C  on 
the  adiabatic  temperature  and  reformate  were  simulated.  From 
Figs.  1-5,  it  can  be  seen  that  the  adiabatic  temperature  and  equi¬ 
librium  compositions  have  a  stronger  dependence  on  the  O/C.  In 
addition,  the  S/C  is  also  a  key  parameter  related  to  the  optimal 
product  compositions. 

A  higher  methanol  conversion  rate  reduces  the  cost,  and 
enhances  the  power  per  volume  or  per  weight  of  the  system.  In 
addition,  any  residual  methanol  can  have  a  harmful  effect  on  the 
proton  exchange  membrane.  Accordingly,  a  complete  methanol 
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Fig.  5.  Effect  of  the  O/C  and  S/C  on  the  mole  fractions  of  C(s)  in  the  reformer  (inlet 
temperature  of  the  reactants  was  573.15  K  at  a  pressure  of  0.1  MPa). 


conversion  is  considered  necessary.  Fortunately,  the  aforemen¬ 
tioned  thermodynamic  analysis  demonstrated  that  methanol  was 
converted  completely  under  the  given  parameters,  as  was  also 
proved  in  [13].  They  additionally  observed  a  100%  methanol  con¬ 
version  rate  at  temperatures  greater  than  623  K. 

The  H2  free- water  molar  fraction  as  a  function  of  the  S/C  and  O/C 
is  depicted  in  Fig.  2.  From  Fig.  2,  it  can  be  observed  that  a  higher  S/C 
favors  H2  production,  and  that  the  H2  effluent  fraction  concentra¬ 
tion  decreases  with  increasing  O/C.  To  exclude  the  impact  induced 
by  the  dilution  of  N2,  a  more  favorable  criterion  should  be  used, 
namely,  the  moles  of  generated  H2  per  mole  of  methanol,  which  is 
also  known  as  the  FI2  yield.  Accordingly,  the  FI2  yield  was  studied 
under  different  operating  conditions,  as  depicted  in  Fig.  3.  When 
the  O/C  was  around  0.32,  a  maximum  value  for  the  FI2  yield  was 
observed. 

To  evaluate  which  intrinsic  factors  caused  equilibrium  compo¬ 
sition  changes  as  a  function  of  the  S/C  and  O/C,  the  equilibrium 
constants  of  all  feasible  reactions  were  plotted  as  a  function  of 
temperature  in  Fig.  8.  It  is  well  known  that,  when  the  equilibrium 
constant,  J<p,  is  much  larger  than  1,  the  reaction  cannot  be  shifted 
by  adjusting  the  reactant  stoichiometric  ratio.  For  a  multi-reaction 
system,  only  the  contribution  of  each  reaction  to  the  product  will 
be  altered;  however  when  I<p  approaches  1,  the  impact  of  reactant 
stoichiometric  ratio  on  the  reaction  widely  diverges.  The  addition 
of  reactants  will  shift  the  equilibrium  to  product  accumulation.  A 


maximized  objective  product  depends  on  the  reactions  for  a  multi¬ 
reaction  system,  whose  I<ps  are  higher.  For  these  reactions,  the 
equilibrium  composition  is  a  strong  function  of  the  temperature. 
For  a  reaction  with  a  I<p  closer  to  1,  the  addition  of  reactants  can 
only  alter  the  relative  magnitude  of  the  maximum.  From  Fig.  1 ,  it 
can  be  observed  that  the  adiabatic  temperature  is  closely  associ¬ 
ated  with  the  O/C,  and  only  a  slight  variation  was  associated  with 
the  S/C.  As  such,  the  maximum  or  minimum  of  the  target  product 
will  primarily  be  influenced  by  the  O/C. 

The  J<pS  associated  to  reactions  involving  H2  are  displayed  as 
a  function  of  temperature  in  Fig.  8.  It  can  be  observed  that  reac¬ 
tions  (1 ),  (2),  (4),  (6)  and  (7)  proceed  more  easily  in  the  investigated 
temperature  range.  Reactions  (10)— (13)  are  influenced  by  equilib¬ 
rium  limitations  in  the  temperature  range  of  700-1000  K.  Reaction 
(5)  is  especially  limited  within  the  entire  investigated  temperature 
range.  At  temperatures  above  1200 K,  reactions  (8)  and  (10)-(13) 
can  be  avoided,  while  their  reverse  reactions  are  more  prone  to 
occur.  Although  reaction  (9)  is  thermodynamically  favorable,  the 
methane  formation  reactions  are  removed,  and  therefore,  reaction 
(9)  will  be  excluded  when  temperatures  exceed  1200  K. 

Reaction  (3)  occurs  much  more  easily  compared  to  other  reac¬ 
tions.  The  Kp  profile  for  reaction  (3)  is  not  depicted  in  Fig.  8,  and 
therefore,  reaction  (3)  is  superior  for  methanol  conversion.  Therein, 
only  excess  methanol  will  be  decomposed  and  steam  reformed  to 
yield  FI2  [14].  Such  reaction  steps  were  proposed  for  the  OSR  of 
methanol  by  Lyubocsky  and  Roychoudhury  [15].  Thus,  a  FI2  fraction 
maximum  exists  when  the  value  of  the  O/C  is  equal  to  0,  as  exhib¬ 
ited  in  Fig.  2;  however,  the  FI2  yield  still  increases  with  increasing 
temperature,  as  depicted  in  Fig.  3.  This  is  the  reason  that  higher  tem¬ 
peratures  suppress  and  eliminate  methanation,  in  reactions  (13) 
and  (12),  and  C  formation,  in  reactions  (11)  and  (10).  The  initial 
reduction  of  FI2  content  with  the  O/C  increase  is  due  to  the  dilution 
of  N2;  however  further  addition  of  air  will  reduce  the  FI2  yield  due 
to  redundant  methanol  combustion.  In  fact,  after  methanation  and 
C  formation  were  efficiently  suppressed  beyond  the  temperature  of 
900  K,  H2  production  is  primarily  ascribed  to  methanol  decompo¬ 
sition  in  reaction  (1),  and  to  methanol  steam  reforming  in  reaction 
(4).  In  comparison  to  the  complete  methanol  oxidation  in  reaction 
(3),  H2  depletion  resulted  from  reaction  (7)  can  be  ignored.  Ideally, 
the  OSR  process  can  be  thermally  neutral  for  an  appropriate  O/C 
and  S/C,  and  the  energy  needed  for  steam  reforming  and  methanol 
decomposition  is  supplied  by  the  complete  methanol  combustion. 
According  to  Fig.  3,  the  favorable  O/C  should  be  between  0.17  and 
0.38  for  a  S/C  of  0-4,  under  which  the  equilibrium  temperature 
should  be  about  1000  K,  as  is  defined  by  area  I  in  Fig.  1. 

C  formation  is  possible  from  reactions  (8)  to  (11).  Although  the 
influence  of  C  formation  on  FI2  yield  can  be  ignored  when  the 
temperature  exceeds  900  K,  these  reactions  are  easily  impacted  by 
operational  parameters,  as  indicated  by  Fig.  8,  due  to  their  rela¬ 
tively  lower  equilibrium  constants.  When  the  temperature  is  lower 
than  1000  K,  a  disproportionation  reaction  is  predominant.  Only 
when  the  temperature  exceeds  1200  K,  are  their  reverse  reactions 
thermodynamically  feasible,  such  that  C  formation  can  be  com¬ 
pletely  removed.  The  corresponding  O/C  is  between  0.4  and  0.58 
for  an  S/C  value  of  0-4,  as  is  referred  to  by  area  II  in  Fig.  1.  In  addi¬ 
tion,  the  S/C  also  changes  the  equilibrium  state  in  the  investigated 
temperature  range  of  700-1200  K.  A  higher  S/C  will  suppress  the 
occurrence  of  reactions  (10)  and  (11).  Therefore,  the  boundary  of 
area  II  is  not  defined  only  as  a  function  of  temperature,  but  is  merely 
a  generalization. 

Fig.  8  also  depicts  the  equilibrium  constants  as  a  function  of  tem¬ 
perature  for  each  reaction  involving  CO.  This  is  to  determine  the 
thermodynamic  role  of  each  reaction  in  CO  formation  and  transfor¬ 
mation.  Reaction  (1)  is  thermodynamically  favorable  throughout 
the  entire  investigated  temperature  range.  Reaction  (3)  is  ther- 
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Fig.  6.  The  equilibrium  composition  and  adiabatic  temperature  profiles  with  respect  to  the  O/C  and  inlet  temperature  at  atmospheric  pressure  when  the  S/C  is  1.5. 


modynamically  the  easiest  occurring  reaction  for  all  reactions 
including  O2,  and  as  a  result,  the  contribution  to  the  CO  reduc¬ 
tion  can  be  neglected  for  reaction  (6).  Reactions  (8),  (10)  and  (12) 
are  favorable  at  low  temperatures,  while  their  reverse  reactions 
are  favorable  at  high  temperatures.  Since  the  I<p  is  always  close 
to  1  within  the  specified  temperature  range  for  reaction  (5),  the 
reaction  is  at  equilibrium  state.  The  addition  of  reactant  can  pro¬ 
mote  the  positive  reaction  or  restrain  the  negative  reaction.  To 


reduce  the  CO  content  as  much  as  possible,  a  higher  S/C  should  be 
selected.  On  the  basis  of  the  above  analysis,  CO  formation  will  be 
enhanced  at  high  temperature;  however  high  temperatures  result 
from  more  methanol  combustion,  which  leads  to  a  reduction  in 
methanol  decomposition.  Thus,  there  exists  a  CO  maximum  with 
increasing  temperature,  as  depicted  in  Fig.  4.  At  the  same  time,  the 
CO  maxima  will  move  to  a  higher  O/C,  accompanied  by  an  increase 
in  the  S/C.  A  CO  maximum  exists  at  an  O/C  of  0.4-0.7,  when  the  S/C  is 
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Fig.  7.  The  equilibrium  composition  and  adiabatic  temperature  profiles  with  respect  to  the  O/C  and  operational  pressure  when  the  S/C  is  1.5  and  the  inlet  temperature  is 
573.15  K. 


between  0  and  4,  and  is  labeled  as  area  III  in  Fig.  1 .  The  correspond¬ 
ing  temperature  is  1300  K,  wherein  the  CO  concentration  reaches  a 
maximum  due  to  the  competition  between  preferential  methanol 
complete  oxidation,  decomposition,  and  the  steam  reforming  of 
excess  methanol. 

An  increase  in  the  S/C  shifts  the  directions  of  reactions  (5) 
and  (10)-(13),  as  they  are  equilibrium  limited.  Therefore,  a  H2 


increase  and  CO  reduction  is  partially  a  result  of  these  reactions 
being  promoted  or  suppressed.  Since  the  occurrence  of  a  H2  and 
CO  maximum  is  primarily  dependant  on  the  temperature,  the  S/C 
can  only  be  of  comparable  magnitude.  In  the  case  of  C  formation, 
this  behavior  is  different.  C  removal  is  simultaneously  dominated 
by  temperature  and  thermodynamic  equilibrium.  When  temper¬ 
ature  exceeds  a  certain  value  or  if  the  S/C  is  sufficiently  high,  C 
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Fig.  8.  The  equilibrium  constants  for  reactions  participating  in  the  methanol  OSR  as 
a  function  of  temperature. 


elimination  will  be  thermodynamically  possible.  This  behavior  is 
also  confirmed  in  references  [4,8]. 


on  the  H2  yield  and  CO  content,  and  therefore,  the  H2  yield  and  CO 
maximum  depend  primarily  on  temperature.  In  addition,  solid  C 
formation  is  derived  from  some  reactions  that  are  reversible  in  the 
investigated  temperature  range.  Hence,  the  elimination  of  solid  C 
is  simultaneously  dominated  by  temperature  and  thermodynamic 
equilibrium.  On  the  basis  of  these  possible  mechanisms,  three  areas 
were  plotted.  Under  area  I,  the  temperature  is  maintained  at  1000  K, 
and  the  relevant  O/C  is  between  0.17  and  0.38  for  a  S/C  of  1-4,  to 
produce  a  maximum  H2  yield.  Under  area  II,  the  boundary  for  solid 
C  formation  is  defined.  When  the  temperature  exceeds  1 200  K,  solid 
C  can  be  removed.  For  area  II,  the  relevant  O/C  is  between  0.4  and 
0.58  for  a  S/C  of  1-4;  however,  this  boundary  still  reverts  to  a  low 
O/C  for  a  higher  S/C  due  to  equilibrium  limitations.  Under  area  III, 
the  CO  content  will  reach  a  maximum,  and  therefore,  the  opera¬ 
tional  temperature  should  be  lower  than  1300  K  to  reduce  the  CO 
molar  fraction  in  the  product  gas.  These  data  provide  guidance  for 
practical  system  development  and  application. 

Acknowledgements 

The  authors  gratefully  acknowledge  the  financial  supports  from 
the  National  High  Technology  Research  and  Development  Program 
of  China  (No.  2007AA05Z148).  We  are  grateful  to  Prof.  S.H.  Chan 
(Nanyang  Technological  University,  Singapore)  for  helpful  advice. 


4.  Conclusions 

A  calculation  method  based  on  FLUENT  software  has  been 
developed  and  successfully  used  to  thermodynamically  analyze 
methanol  autothermal  reforming  for  PEMFC.  The  effects  of  temper¬ 
ature,  pressure,  S/C,  and  O/C  on  the  OSR  of  methanol  were  identical 
to  previous  published  results.  Slight  variations  in  the  data  are  due  to 
differences  in  the  thermodynamic  database  used  in  the  two  meth¬ 
ods.  The  equilibrium  constants  and  the  possible  reactions  involved 
in  OSR  of  methanol  were  calculated  to  elucidate  the  reaction  mech¬ 
anisms  in  the  entire  investigated  temperature  range.  It  can  be 
concluded  that  four  major  reactions,  including  methanol  combus¬ 
tion,  methanol  decomposition,  steam  reforming,  and  water-gas 
shift  reaction,  are  involved  in  OSR  of  methanol.  Methanol  com¬ 
bustion  was  observed  to  be  dominant,  while  only  excess  methanol 
decomposed  and  steam  reformed.  Reversible  water-gas  shift  reac¬ 
tions  also  occur  in  the  investigated  temperature  range.  Further, 
other  reversible  reactions  involving  H2  and  CO  have  a  slight  effect 
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